Manipulation of the Lande g-factor in InAs quantum dots through the application of 
anisotropic gate potentials: Exact diagonalization, numerical and perturbation 

methods 
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We study the variation in the Lande g-factor of electron spins induced by both anisotropic gate 
potentials and magnetic fields in InAs quantum dots for possible implementation towards solid state 
quantum computing. In this paper, we present analytical expressions and numerical simulations of 
the variation in the Lande g-factor for both isotropic and anisotropic quantum dots. Using both 
analytical techniques and numerical simulations, we show that the Rashba spin-orbit coupling has a 
major contribution in the variation of the g-factor with electric fields before the regime, where level 
crossing or anticrossing occurs. In particular, the electric field tunability is shown to cover a wide 
range of g-factor through strong Rashba spin-orbit interaction. Another major result of this paper 
is that the anisotropic gate potential gives quenching effect in the orbital angular momentum that 
reduces the variation in the E-field and B-field tunability of the g-factor if the area of the symmetric 
and asymmetric quantum dots is held constant. We identify level crossings and anticrossings of 
the electron states in the variation of the Lande g-factor. We model the wavefunctions of electron 
spins and estimate the size of the anticrossing for the spin states |0, —1, +1/2 > and |0, 0, —1/2 > 
corresponding to a quantum dot that has been recently studied experimentally (Phys. Rev. Lett. 
104, 246801 (2010)). 
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I. INTRODUCTION 

Single electron spins in an electrostatically defined 
quantum dot in a 2-dimensional electron gas (2DEG) 
have been manipulated and studied by several groups 
Quantum dots in III-V type semiconductors provide an 
opportunity to study the variation in the Lande g-factor 
vs. gate potentials and magnetic fields^— The shape 
and size of the quantum dots can be modified by chang- 
ing the gate controlled electric fields that influence the 
variation in the energy spectrum as well as the Lande 
g-factor of the dots^ ' 14 i 15 The results of this research 
might enhance the opportunities of building spintronic 
logic devices for possible implementation towards solid 
state quantum computing 

The orbital and spin angular momentum of the elec- 
tron in a semiconductor quantum dot interact through 
the Rashba and Dresselhaus spin-orbit couplingSi 21 i 22 
These two spin-orbit coupling effects arise from two dif- 
ferent types of symmetry operations in III-V type semi- 
conductors. The Rashba spin-orbit coupling arises from 
the structural inversion asymmetry of the triangular 
shaped quantum well confining potential. The Dressel- 
haus spin-orbit coupling arises from bulk inversion asym- 
metry. The mathematical expressions for the Rashba and 
Dresselhaus spin-orbit couplings that are implemented 
into the theoretical model are given in this paper in 
Section II. The strength of the Rashba and Dresselhaus 
spin-orbit couplings is determined by the gate controlled 
electric fields and is an important parameter in control- 
ling the Lande g-factor for both isotropic and anisotropic 
quantum dots. 

The goal of the present work is to explore the non- 



degenerate energy spectrum of electrostatically defined 
InAs quantum dot by both analytical techniques and 
numerical simulations. These approaches provide real- 
istic information for controlling the Lande g-factor for 
both isotropic and anisotropic quantum dots through the 
application of gate potentials. We also model realistic 
wavefunctions of electrons in InAs quantum dots that 
were recently studied by experimentalists in Ref. I23I . 
We estimate the size of avoided anticrossing is approx- 
imately 65 /xeV, which is in agrement with the exper- 
imentally reported values* 23 In this paper, by utilizing 
both analytical and numerical techniques, we find that 
the Rashba spin-orbit coupling produces the dominant 
effect on the variation of the Lande g-factor vs. elec- 
tric field strength below the level crossing or anticross- 
ing. Also, anisotropic gate potential lead to a quenching 
effect in the orbital angular momentum that reduces the 
variation in the g-fa ctor. Our work is similar to those of 
Refs. S, 0, H, [H IH and [H but differs in that we utilize 
both analytical and numerical approaches based on the 
finite element method. 

The paper is organized as follows. In Sec. II, we ex- 
actly diagonalize the Hamiltonian of a quantum dot con- 
fined in an asymmetric potential including spin-orbit in- 
teraction and a magnetic field along z-direction. In Sec. 

III, based on a second order perturbation calculation, 
we present analytical expressions of the Lande g-factor 
for both isotropic and anisotropic quantum dots. In Sec. 

IV, we plot the Lande g-factor induced by gate controlled 
electric fields vs. magnetic field as well as quantum dot 
radii for pure Rashba case (ccd = 0), pure Dresselhaus 
case {an — 0) and mixed cases (both an,ao present) 
where an and old are two parameters related to the 
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II. THEORETICAL MODEL 

The Hamiltonian of an electron in a quantum dot in the 
plane of a 2DEG, in the presence of an external magnetic 
field along z-direction can be written as^l 



H — H r 



(1) 



where the Hamiltonian H so is associated to the Rashba- 
Dresselhaus spin-orbit couplings and H xy is the Hamil- 
tonian of the electron along the lateral direction in the 
plane of the 2DEG. H xy can be written as 



p2 2 

H xy = — + -mu 2 (ax 2 



by 2 



-^9oliBO z B, (2) 



FIG. 1. Contributions of Rashba and Dresselhaus spin-orbit 
couplings {cur/ old) vs - Electric Field. It can be seen that 
Rashba and Dresselhaus spin-orbit couplings becomes equal 
at the electric field, E = 3.05 x 10 3 V/cm. This plot is ob- 
tained from Eq. |9] 



strength of the Rashba and Dresselhaus spin-orbit in- 
teractions respectively. We also model the wavefunctions 
of electron states that were recently reported by experi- 
mentalists in Ref. l23l . Finally, in Sec. V, we summarize 
our results. 
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where P = p + eA is the kinetic momentum operator, 
p = —ih(d x ,d y ,0) is canonical momentum operator and 
A = -M^{—yb, xa, 0) is the vector potential in the asym- 
metric gauge. Here, — e < is the electronic charge, m is 
the effective mass of the electron in the conduction band, 
fiB is the Bohr magneton, a z is the Pauli spin matrix 
along z-direction. Also, uiq — is a parameter charac- 
terizing the strength of the confining potential and £q is 
used as the radius of the quantum dot. By choosing the 
appropriate values of a and 6, one can define the shape 
of the quantum dot from circle to ellipse. The Hamil- 
tonian H xy can be exactly diagonalized (see Appendix 
A) by writing the Hamiltonian as a function of number 
operators n± = a±a± and its energy spectrum can be 
written asi2£ 



a± 



1 



(«+-*_) + 
1 

(*+-0(l-i) 



£ n+n ^ = (n+ + n- + 1) huj + + (n + - n_) fwj- + -goHB<?zB, 

I I 1 1 " 

± (s T ± i) tPx + (s± ± i) -p>y + (1 =F is T ) jX±(It is±) jV 

I i 1 1 

± (s T =F i) jVx + (s± T i) -zPy + (1 ± is T ) jx±(l± is±) -y 



(3) 
(4) 

(5) 



where, 




" ±= 2 



1/2 



provided that 



a±,a± 



1 and [x,p x ] = [y,p y ] 



(6) 

(7) 
it). 



Here fij = uj ay /a, fi 2 = u Vb, £ = J ^ 



+ and uj c — — is the cyclotron frequency. 
Note that a similar type of expression for the energy 
spectrum of an anisotropic quantum dot is also discussed 



in Ref. |26| . However, our methodology for diagonaliz- 
ing the Hamiltonian is slightly different (see in Appendix 
A). Also, we verified that by substituting a — b — 1, Eqs. 

are exactly the same as those in Ref. 0- The eigen- 
states of H xy in Eq. [5] with a = b = 1 arc well-known 
Fock-Darwin energy states . 29 ' 30 

The Hamiltonian associated with the spin-orbit cou- 
plings can be written a a 21 i 22 

H so = -T- \PxPy - VyPx) + ~J~ (-CT X P X + CTyPy) ■ (8) 

The spin-orbit Hamiltonian consists of the Rashba cou- 
pling whose strength is characterized by the parameter 
an and the Dresselhaus coupling with olu- These cou- 
pling parameters depend on the electric field E of the 
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quantum well confining potential (i.e., E = —dV/dz) 
along z direction at the interface in a hetero junction as 



/ <2yjk> \ 2 ^ 

a R = lR eE, a D =0.78 lD (^J £ 2/3 , 



(9) 



J 



where the Rashba coefficient 7^ = llOA , Dresselhaus 

o 3 

coefficient 7^ = 130 eVA and effective mass m = 0.0239 
have been considered for InAs quantum dots^ 

The Rashba and Dresselhaus spin-orbit Hamiltonian 
can be written in terms of raising and lowering operators 
as>i 



H so = a R (1 + i) b x ^ 4 n + (s_|_ — i) a + + b x ^ 4 n + (s_ + i) a_ + a 1 / 4 ^- (i — s_) a + + a 1 / 4 ?]- (i + ,s + ) a_ 

+£*£>(! + i) a 1/4 K- (i - s-)a+ + a 1/4 K^ (i + s+)a- +b 1/4 r] + (-i + s+)a+ +b 1/4 r]+ (i + s-)a- + h.c, (10) 



where. 



1 [1 , . eB£ 

; r s —<J r ± 1 

2(s+-s-)\e n 

1 Jl .eS^ 

2(s+- S _) jl^ ~?T 
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v/a+ \/6 
1 



'2/ r > 



and h.c. represents the hermitian conjugate. It is clear 
that the spin-orbit Hamiltonian and the Zeeman spin 
splitting energy in both isotropic and anisotropic quan- 
tum dots obeys a selection rule in which the orbital an- 
gular momentum can change by one quantum. 



III. ANALYTICAL EXPRESSION FOR THE 
LANDE j-FACTOR 

Our conventional definition of the electron g-factor in 
the conduction band in the presence of magnetic field 
along z-direction can be written as^i 



e 0,0, + l/2 — £ 0, 0,-1/2 

fi B B 



(11) 



where £0,0, +1/2 an d £ o, 0,-1/2 are the eigenvalues of the 
electron in the conduction band with spin up and down 
respectively having the lowest orbital angular momen- 
tum. 

At low electric fields and small quantum dot radii, we 
treat the Hamiltonian associated to the Rashba and Dres- 
selhaus spin-orbit couplings as a perturbation. Based on 
second order perturbation theory, the energy of the elec- 



r 



tron spin states can be written as, 

_ (0) , (2) 

£ 0,0, + l/2 — £ 0,0, + l/2 + e : 0, + l/2' 

_ (0) , (2) 

£ 0,0,-l/2 — £ 0, 0,-1/2 + £ 0,0,-l/2' 



(12) 
(13) 



The zero order energy correction can be easily calculated 
from Eq. ^ . The first order energy correction is zero. 
The calculation for the second order energy correction 
can be written as: 



e {2) 

c 0,0,+l/2 

(2) 
e 0, 0,-1/2 

where, 

£± = 
c± = 



a 2 R £+ + a 2 D e; + 



jo) 

-0,0, +1/2 



-(0) 

'1,0,-1/2 



(0) _ (0) 

£ 0, 0,-1/2 £ 1, 0,1/2 



JO) _ JO) 

t 0,0, + l/2 fc 0, 1,-1/2 



a|g_ + a 2 D <;- 
-(0) _ (o) 

'0,0,-1/2 £ 0,l, + l/2 



,(14) 



,(15) 



1 



2( S+ - S _) 
1 

2(s+ - s-) 



±— al+2a ± f3 ± T— 0± ^,(16) 

S± Sz, 

1 2 
s± 



2a T P T T —Pi ,(17) 



1/4 J 1 1 eB£ 
a± = a ' < - ± 



h (V^+Vb) 



(y/E+Vb 



,(18) 



(19) 



By substituting Eqs. [T2l and [TBI in Eq. QTJ the expression 
for the Lande g-factor of anisotropic quantum dots can 
be written as, 



J 



b i^bB \ huj x — A hujy — A huj x + A hujy + A ' 



I 

In the above expression, we use the relation u> x = u> + + oj_, uj y = u> + — u;_ and A = q^^bB. By substituting 
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b = 1 in Eq. 1201 one find the Lande g- factor for isotropic quantum dot: 
I 
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[al(l-8)-a 2 R (l + 5)] I 2 
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(21) 
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FIG. 2. (Color online) (a) Eigenenergy levels for ground, first 
and second excited states vs. magnetic fields, (b) Effective 
Lande g-factor vs. magnetic fields. In both cases, solid lines 
represent pure Dresselhaus case (as = 0), dashed lines rep- 
resent pure Rashba case (an — 0) and open circles represent 
mixed cases (both olr and old present). Also in both cases, 
we choose the electric field E = 3 x 10 5 V/cm and the quan- 
tum dot radius, £q = 20 nm for symmetric quantum dots 
(a = 6 = 1). 



where go — —15 is the bulk g-factor for InAs quantum 
dots and 5 — gQm/m e . Here m e is the mass of electron. 
The analytical expressions for the Lande g-factor of both 
anisotropic and isotropic quantum dots in Eqs. (|20p and 
(f2"Tj) respectively are valid before the level crossing or an- 
ticrossing occurs. It can be seen that the g-factor for 
both isotropic and anisotropic quantum dots depends on 
the anisotropic gate potentials, quantum dot radii and 
magnetic fields. These are our control parameters in ma- 
nipulating the g-factor for both isotropic and anisotropic 
quantum dots. 



IV. RESULTS AND DISCUSSIONS 

We have used the Finite Element Method 31 to solve 
numerically the corresponding eigenvalue problem with 
Hamiltonian given by Eq. ([T]) to study the variation in 
the Lande g-factor at high electric fields for the values 
olr/ old > 1, large quantum dot radii and magnetic fields 
for both isotropic (a = b) and anisotropic (a ^ b) quan- 
tum dots. Throughout the simulations, unless otherwise 
stated, we consider £q — 20 nm for g/gQ vs. magnetic 
field and B = 1 T for g/go vs. quantum dot radius. Our 
numerical simulations related to the g-factor in quantum 
dots are valid before the level crossing or anticrossing 
occurs. At or after the level crossing or anticrossing, 
g-factor simply captures the effects of two lowest eigen- 
states in quantum dots. 

In Fig. [TJ we investigate the strength of the Rashba 
and Dresselhaus spin-orbit interactions vs. electric fields. 
The strength of the Rashba and Dresselhaus spin-orbit 
couplings for InAs quantum dots can be determined by 
the relation ajj/ac = 0.015-E 1 / 3 (see Eq. [5]). From this 
relation, it can be found that both Rashba and Dres- 
selhaus spin-orbit couplings become equal at the electric 
field E = 3.05 x 10 3 V/cm. However, at this value of the 
electric field where olr = ajj, the rotational symmetry 
is not broken and the spin splitting energy mainly corre- 
sponds to the Zeeman energy. For E = 10 4 V/cm to 10 6 
V/cm, an/ old varies from 1.5 to 6.88 (see Fig.QJ. Since 
olr/old > 1, only the Rashba spin-orbit coupling has an 
appreciable contribution to spin splitting energy. 

Figure [2] illustrates the eigenenergy of the states 
|0,0,+l/2 >, |0, 0,-1/2 >, and |0,-l,+l/2 > and g- 
factor vs. magnetic fields of symmetric quantum dots 
(a = b = 1) for pure Rashba {old = 0, solid lines), 
pure Dresselhaus (olr = 0, dashed lines) and mixed cases 
(both olr and old present, open circles). It can be seen 
that the Dresselhaus spin-orbit coupling has almost no 
effect on the manipulation of the g-factor. For the pure 
Dresselhaus case, we find a level crossing at 3.5 T. How- 
ever, for the pure Rashba and mixed spin-orbit couplings, 
we find an avoided anticrossings in the manipulation of 
the g-factor. 

Figure [3] explores the variation of several eigenenergy 
states and the g-factor vs. magnetic field for a symmetric 
quantum dot (a = b = 1). In Fig. [3]Ja), we plot several 
eigenenergy levels (ground, first, etc.) vs. magnetic field 
for symmetric quantum dots with gate induced electric 
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fields E = 10 4 V/cm (solid lines) and E = 5 x 10 5 V/cm 
(dashed lines). Here we find level crossing shown by open 
circles at magnetic field, B = 3.5 T for the electric field 
E = 10 4 V/cm (solid lines) . However, we find an anti- 
crossing for electric field E = 5xl0 5 V/ cm (dashed lines) . 
The symbol x in Fig. [3] (a) represents the data from per- 
turbation theory which is in agreement with the numeri- 
cal simulations. In Fig. [2Kb), we plot the Lande g-factor 
vs. magnetic field at the electric fields E = 10 4 V/cm 
(solid lines), 3 x 10 5 V/cm (dashed lines), 5 x 10 5 V/cm 
(dotted lines), 7 x 10 5 V/cm (dashed-dotted lines) and 
10 6 V/cm (dashed-dotted-dotted lines) for a symmetric 
quantum dot. It is clear that there is a level crossing be- 
tween the states |0, 0, -1/2 > and |0, -1, +1/2 > in the 
manipulation of the g-factor at B = 3.5 T for E = 10 4 
V/cm (solid lines). However, with the increase in the 
electric field, we have an anticrossing in the variation of 
the g-factor. This is also consistent with previously pub- 
lished work by R. de Sousa and S. Das Sarma£ and we 
consider this result as a benchmark for our computational 
work. 

In Fig. HI we analyze anisotropy effects on the variation 
of the eigenenergy and of the Lande g-factor vs. magnetic 
fields. In Fig.QJa), we plot the eigenenergy vs. magnetic 
field for quantum dots in the potentials characterized by 
a = b = 1 (solid lines), a = b = 3 (dashed lines) and a = 
1,6 = 9 (dotted lines). Here, we choose E = 10 4 V/cm. 
We find a level crossing at B — 3.5 T for the symmetric 
quantum dot (a = b = 1). However, the level crossing 
extends approximately to B = 6.2 T for the symmetric 
quantum dots (a = b = 3). The extension of the level 
crossing to the larger magnetic fields is mainly due to 
the increase on the lateral size of the quantum dots. We 
quantify that the area of the quantum dots with (a = b = 
3) is 9 times larger than the dots with (a = b = 1). For 
anisotropic quantum dots with the potential (a = 1,6 = 
9), we find the level crossing at around B — 6 T. Note 
that the area of the asymmetric dots in the potential (a = 
1 , 6 = 9) is exactly equal to the area of the symmetric dots 
with the potential (a = 6 = 3). In Fig. H](b) and (c), we 
plot the variation in the Lande g-factor vs. magnetic field 
for quantum dots in the confining potentials of a = b = 3 
and a = 1,6 = 9. In both cases, we choose the electric 
fields as same as in Fig. [3Kb). In Fig.@](b) for symmetric 
quantum dots (a=b=3), we find a crossover point (all 
the curves collapse to a single point) at around B = 5 T. 
However, this point decreases to approximately B = 4.5 
T for asymmetric quantum dots (a=l, b=9) as shown in 
Fig. Hfc). At this point, the value of g-factor is equal to 
the bulk g- value of InAs dots (i.e., g = go = — 15) for 
all gate induced electric fields. Also, in Fig. [He), we see 
that the anisotropic potential gives a quenching effect in 
the orbital angular momentum that reduces the electric 
field tunability of the g-factor. 

Finally, Fig. [5] illustrates the manipulation of the 
Lande g-factor vs. quantum dot radii for both symmet- 
ric and asymmetric quantum dots. Again, the electric 
fields are chosen as same as in Fig.[3Kb). In Fig.fSJa) and 
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FIG. 3. (Color online) (a) Several eigenenergy levels vs. mag- 
netic fields at the electric fields: E = 10 4 V/cm (solid lines) 
and E = 5 x 10 V/cm (dashed lines). The symbol x repre- 
sents data from perturbation theory. A level crossing occurs 
at the magnetic field, B = 3.5 Tesla for the electric field, 
E = 10 4 V/cm (solid line) and is shown by open circles, (b) 
Effective Lande g-factor vs. magnetic field at the electric fields 
E = 10 4 V/cm (solid line), 3xl0 5 V/cm (dashed line), 5xl0 5 
V/cm (dotted line), 7 x 10 5 V/cm (dashed-dotted line) and 
10 6 V/cm (dashed-dotted-dotted line). In (a) and (b), we 
choose the quantum dot radius, £q — 20 nm for symmetric 
quantum dots (a = b = 1). 



(b), we consider isotropic confining potentials: a = b = 1 
and a = b = 3. In Fig. [5th), level crossing point ex- 
tends to the larger quantum dot radii compared to that 
of Fig. EJa) mainly due to the increase in the lateral 
size of the quantum dots. In Fig. [SJc), we capture the 
anisotropic effect of the E-field tunability of the g-factor 
vs. quantum dot radii with the potential a = 1,6 = 9. 
Again, we see that the anisotropic potential a = 1,6 = 9 
(Fig. (5](c)) reduces the E-field tunability compared to 
that of isotropic potential a = b = 3. 

Now we consider the wavefunctions and estimate the 
size of the anticrossing point in an experimentally re- 
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FIG. 4. (Color online) (a) Several eigenenergy levels vs. magnetic field at the electric field E = 10 4 V/cm for both symmetric 
(a — b — 1 (solid lines), a — b = 3 (dashed lines)) and asymmetric quantum dots (a = 1, b — 9, (dotted lines)). In (b) and (c), 
we plot effective Lande g-factor vs. magnetic fields. The electric fields are chosen as same as in Fig. O^b). Also, in (b) and (c), 
we choose a — b = 3 and a = 1, b = 9 respectively. In all cases, we choose the quantum dot radius, £q = 20 nm. 
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FIG. 5. (color online) Lande g-factor vs. quantum dot radius for both symmetric (a = b) and asymmetric (a 7^ b) quantum 
dots. Again, the electric fields are chosen as same as in Fig. [3jb). Also we choose = 6= 1, = 6 = 3 and a = 1, 6 = 9 in 
(a),(b) and (c) respectively. In all cases, we consider the magnetic field B = 1 T. 
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FIG. 6. (color online) In-plane wavefunctions in an asymmetric quantum dot. Here we choose a = 1.5, b — 4, 
E = 1.6 X 10 4 V/cm, B = 2.9 T and £ = 28 nm. These parameters mimic the wavefunctions of the quantum states 
|0,0,+l/2 > , |0, -1.+1/2 > and |0,0, -1/2 > (from left to right) in Ref. [H of Fig. (2a). 
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FIG. 7. Eigenenergy of electron states |0, 0, +1/2 > 
, |0, —1, +1/2 > and |0, 0,-1/2 > vs. magnetic field. The 
symbol x represents the data from perturbation theory. The 
inset shows the energy difference vs. magnetic field near the 
level crossing point. Here we estimate the size of the anti- 
crossing point to be approximately 65 \ieV . Here we again 
consider a = 1.5, 6 = 4, E = 1.6 X 10 4 V/cm, and £ = 28 
nm. 
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FIG. 8. (Color online) Manipulation of the Lande g-factor vs. 
magnetic and electric fields for quantum dots in the potential 
characterized by (top to bottom) a — b = 3 (solid lines), a = 
2, b = 4.5 (dashed lines), a — 1.5,6 = 6 (dashed-dotted line) 
and a = 1,6 = 9 (dashed-dotted-dotted lines). We consider 
electric field, E = 7 x 10 5 V/cm in case (a) and magnetic 
fields B = 1 T in case (b). In both cases, we consider the 
quantum dot radius to = 20 nm. 



ported self assembled InAs quantum dot in Ref. [23|. 
The authors in Ref. [23| characterize the lateral size of 
the quantum dots by two anisotropic gate potentials, 
huj x — 1.5 meV and hcu y — 4 meV. This anisotropic 
potentials give 46 nm lateral size of the quantum dots 
along the x-axis and 28 nm along the y-axis in the plane 
of 2DEG. By substituting a = 1.5, 6 = 4 and £ = 28 nm 
in our theoretical model, this will mimic the experimen- 
tally reported values of the lateral size of the quantum 
dots. The authors in Ref. [H also report that the height 
of the pyramidal shape quantum dots is 20 nm. In our 
theoretical model, the quantum dot is formed in the plane 
of 2DEG so we choose the average height of the quantum 
dots as 10 nm. In very crude approximation, this thick- 
ness of the 2DEG is their height of the pyramidal shape 
quantum dots^ Now, the expression (2meE /?i 2 ) _1 / 3 in 
our theoretical model gives an estimate of the thickness 
of 2DEG which basically defines the size of the quantum 
dot along z-direction. 10 nm vertical size of the quan- 
tum dot in the above expression gives an estimate of the 
electric fields, E = 1.6 x 10 4 V/cm. We find the anti- 
crossing point at around B = 2.9 Tesla and consider this 
value to capture the realistic wavefunctions for the states 
|0, 0, +1/2 >,|0, 0,-1/2 > and |0, -1,1/2 >. From Fig. 
the first excited state wavefunction is |0, — 1, +1/2 > 
which is a clear indication of the level crossing point. 
In Fig. [7J we plot eigenenergy of the electron states 
|0, 0+1/2 >, |0, —1, +1/2 > and |0, 0,-1/2 > vs. mag- 
netic field. Again, the symbol x in Fig. [7J represents the 
data from perturbation theory for asymmetric potentials 
which is in agreement with the numerical simulations. 
It can be seen that there is an anticrossing point at the 
magnetic field of around 2.9 Tesla. Theoretically inves- 
tigated level crossing point (B=2.9 T) is slightly smaller 
than the experimentally reported values (B=3.5 T, see 
Ref. l23f ) which indicates that the characterization of the 
lateral size of the quantum dots by anisotropic potentials 
has an estimate error of 20%. The inset plot shows the 
magnified image near the level crossing point. Here we 
estimate the size of the anticrossing to be nearly 65 /zeV 
which is in a gree ment with the experimentally reported 
values in Ref. [23|. 

In Fig. [8l we summarize our results for the manipu- 
lation of the Lande g-factor with respect to both mag- 
netic and electric fields in isotropic and anisotropic quan- 
tum dots. In Fig. Ufa), we plot the g-factor vs. mag- 
netic field at E = 7 x 10 5 V/cm. Also, we plot the 
g-factor vs. electric field at B = 1 T in Fig. \%lb). In 
both cases, we compare the result for isotropic quantum 
dots (a = b = 3, solid lines) and for anisotropic quan- 
tum dots a = 2,b = 4.5 (dashed lines), a = 1.5,6 = 6 
(dashed-dotted lines) and a = 1,6 = 9 (dashed-dotted- 
dotted lines). Note that the quantum dots in the above 
confining potentials have the same area. At large mag- 
netic fields, B ~ 5 T in Fig.[8ja) and small electric fields, 
E ~ 10 5 V/cm in Fig. EKb) , the effect of the resulting spin 
and orbital angular momentum of the quantum states on 
the eigenvalues (eo,o,+i/2 — £0,0,-1/2) are small, leading 
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towards the bulk g-factor i.e. g/go ~ 1 or g = —15. How- 
ever, at small magnetic fields, B ~ 1 T in Fig. [SJa) and 
large electric fields, E ~ 10 6 V/cm in Fig.[8jb), the effect 
on (e ,o,+i/2 - £0,0,-1/2) is large enough to drag the g- 
factor towards the free electron value i.e. g/go ~ 2. Also, 
in both cases, if we compare the g-factor in both isotropic 
and anisotropic quantum dots, we find that anisotropic 
potentials (a + 1 b) give the quenching effect in the or- 
bital angular momentum that reduces the variation in 
the B-field tunability and E-field tunability of the Lande 
g-factor. This is also expected from Eqs. I2U1 and |2"T1 For 
example, suppose B — 1 T, £ — 20 nm and E = 3 x 10 5 
V/cm. Then from Eq. [20j we find g aS ym = 1-04 for 
anisotropic quantum dots (a = 1, b = 9). However, from 
Eq. 1211 we find g S ym = 1-07 for isotropic quantum dots 
( a = b = 3). Also in Fig. [S] (a) and (b), it can be seen 
that g sy m > gasym- Again, this quantitative analysis is 
in excellent agreement with our numerical simulations. 



V. CONCLUSIONS 

By utilizing both analytical expressions and numeri- 
cal simulations (based on Finite Element Method), from 
Figs.Q~|[51 we have shown that the Rashba spin-orbit cou- 
pling has the major contribution on the variation of the 
g-factor with electric fields before the level crossing or 
anticrossing occurs. In Fig. [6l we have demonstrated 
the results of modeling of the realistic wavefunctions of 
the states |0, 0, +1/2 >, |0, — 1, +1/2 > and |0, 0,-1/2 > 
near the level crossing point and estimate 65 ^eV as the 
size of the anticrossing in Fig. [7] Finally, in Fig. |8l we 
have shown that the anisotropic gate potential gives the 
quenching effect in the orbital angular momentum that 
reduces the variation in the B-field tunability and E-field 
tunability of the g-factor. 
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X 2 



V, 



P> 



^0 

fi 2 y ' 



(A3) 



where f^i = wov^ an d ^2 = ^oVb. Also, the Gauge 
potential can be written as 



A, 



x 2 B^/il 2 uj 

n~i + n 2 ' 



a,< 



XiBy/TlxLOQ 



(A4) 



Substituting Eqs. (JAUEHIEII into Eq. (|AT) . we get the 
Hamiltonian in the form of: 



[Pi + x\ + e (Pf + xj) + c { Xl P 2 - x 2 P 1 )] 



(A5) 

5| are as follows: 

1/2 



lo 2 c + (fil + Q 2 f 



2 m wo 

The abbreviations used in Eq. 
e = fi 2 /ni, c = (2w cV /0 2 /^i) / 
Also we use the relation mcjo7 = 1, where = 1 + 

The energy spectrum of Hamiltonian (|A5|) can be 
found as follows. First, we need to find the canonical 
transformation U of the four-dimensional phase space, 
P = (P x , P y , x, y) which diagonalizes the quadratic form 
of the Hamiltonian (|A5[) . To be more specific, Hamilto- 
nian (lAll) can be written as 



-c/2 
e c/2 



H - = (2^) PtMP ' M 




(A6) 

where t represents the transpose of a vector. The or- 
thogonal unitary matrix U which exactly diagonalizes the 
matrix M can be written as, 



/ 1 



U 



(*+-*-) 



1 -1 

S- 



— S- 


-s+ 


s+ 


— S- 


1 


1 


1 


-1 



(A7) 
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Appendix A: Energy spectrum of asymmetric 
quantum dot 

The Hamiltonian of an electron in a quantum dot in 
the plane of a 2DEG can be written as 



<-±^L + l - m ul(ax 2 +by 2 ). (Al) 



The Eq. (|A1[) can be exactly diagonalized^ by making 
a canonical transformation of position and momentum 
operators as 



Xi 




Pi 



(A2) 



1 ± d and 




(A8) 



Also, the expression for s± is written in Eq. (j6|). In terms 
of rotated operators P' = UP, the Hamiltonian (IA6|) can 
be written as 



H , 



2tocl>o 



-( CS _+2)(pl 2 



y 



'2 



The above Hamiltonian is identified as the superposition 
of two independent harmonic oscillators and its energy 
spectrum can be written as 

e n+ n_ = (n+ + n_ + 1) hcu + + (n+ - n_) ftcj_, (A9) 
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provided that 



a±,a± 



= 1. 



Also ui± = 



1/2 
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